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A  MODEL  OF  THE  EFFECTS  OF  PRESSURE  AND  CROSSFLOW  VELOCITY  ON 
COMPOSITE  PROPELLANT  BURNING  RATE* 

Dr.  Merrill  K.  King** 
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Abstract 


Several  variations  of  a  model  for  prediction  of 
burning  rate  versus  pressure  behavior  of  unimodal 
oxidizer  composite  propellants  in  the  absence  of 
crossflow  were  developed  and  evaluated  against  a 
set  of  data  for  a  series  of  four  formulations. 

Three  variants,  including  one  in  which  an  average 
oxidizer-burning-surface  intersectional  area  con¬ 
cept  is  employed  and  two  in  which  allowance  is  made 
for  geometry  and  stoichiometry  changes  as  the  pro¬ 
pellant  recedes  past  an  oxidizer  crystal,  were 
found  to  give  excellent  agreement  with  data.  The 
former- variant  was  ejftended  to  treat  multimodal 
oxidizer  formulations,  yielding  predictions  in 
excellent  agreement  with  data  for  two  additional 
formulations  containing  bimodal  oxidizer.  In  the 
Initial  development  of  the  erosive  burning  aspect 
of  the  model,  columnar  diffusion  flame  bending  was 


assumed  to  be  the  sole  mechanism  leading  to  burning 
rate  augmentation  by  crossflow:  this  assumption 
led  to  severe  underprediction  of  erosive  burning 
effects.  Accordingly,  the  model  was  revised  through 
addition  of  a  flow  profile  analysis  for  prediction 
of  cross-flow-induced  turbulence  augmentation  of 
transport  properties  governing  heat  feedback  from 
gas  flames  as  well  as  flame-bending.  This  revised 
model  was  found  to  yield  good  agreement  with  ero¬ 
sive  burning  data  for  five  of  the  six  formulations 
tested,  but  gave  higher  predicted  rates  than  ob¬ 
served  values  for  the  sixth  propellant. 

Introduction  and  Background 

Erosive  burning,  the  augmentation  of  solid  pro¬ 
pellant  burning  rate  by  the  flow  of  products  across 
a  burning  surface,  is  becoming  increasingly  import¬ 
ant  with  use  of  lower  port-to-throat  area  ratio 
motors  and  nozzleless  motors,  both  of  which  result 
in  high  velocity  crossflows.  The  response  of  var¬ 
ious  propellants  to  such  crossflows  must  be  known 
by  the  motor  designer  in  order  for  him  to  perform 
adequate  motor  design.  In  addition,  it  is  import¬ 
ant  that  the  propellant  formulator  understand  the 
effect  of  various  formulation  parameters  on  the 
sensitivity  of  a  propellant  to  crossflows  so  that 
he  nay  tailor  his  propellants  to  the  desired  char¬ 
acteristics.  For  example,  in  a  nozzleless  rocket 
motor,  the  decrease  in  pressure  from  the  head  end 
to  the  aft  end  of  the  grain  tends  to  result  in 
slower  burning  at  the  aft  end  in  the  absence  of 
erosive  effects.  Depending  upon  the  sensitivity 
of  the  formulation  to  crossflow,  the  increasing 
Mach  Number  along  the  grain  port  may  lead  to  under¬ 
compensation,  exact  cancellation,  or  over-compensa¬ 
tion  of  the  pressure  effect.  A  detailed  discussion 
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of  the  effects  of  erosive  burning  on  solid  propel¬ 
lant  rocket  interior  ballistics  for  low  port-to- 
throat  area  ratio  motors  and  nozzleless  motors  was 
presented  by  this  author  in  Reference  1.  During  the 
past  three  years,  this  author  has  been  conducting 
an  analytical  and  experimental  study  of  the  erosive 
burning  of  composite  propellants.  This  effort  in¬ 
cludes:  (1)  development  of  a  simplified  model  for 
prediction  of  the  erosive  burning  of  a  composite 
propellant,  given  the  non-erosive  burning  rate  ver¬ 
sus  pressure  relationship  for  that  formulation; 

(2)  development  of  a  more  fundamental  composite  pro¬ 
pellant  combustion  model  for  prediction  of  burning 
rate  as  a  function  of  pressure  and  crossflow  vel¬ 
ocity,  given  only  the  propellant  composition  and  in¬ 
gredient  particle  size  distributions;  (3)  experi¬ 
mental  measurement  of  the  erosive  burning  charact¬ 
eristics  (at  crossflow  velocities  up  to  Mach  1)  of  a 
series  of  propellants  with  systematically  varied 
compositions  and  Ingredient  particle  sizes;  and  (4) 
fine  tuning  of  the  models  using  this  experimental 
data.  The  simplified  (Generation  1)  model  has  been 
described  in  detail  in  References  2  and  3,  with  the 
experimental  procedure  and  test  results  to  date 
appearing  in  References  3  and  4. 

Details  of  the  first  variant  of  the  more  funda¬ 
mental  (Generation  2)  model  for  the  prediction  of 
burning  rate  as  a  function  of  pressure  and  crossflow 
velocity  appear  in  Reference  5.  (That  reference 
will  be  cited  heavily  throughout  this  paper,  with 
later  variants  of  the  Generation  2  model  being  des-  • 
cribed  in  terms  of  modifications  of  the  analysis 
presented  therein.)  This  first  variant  of  the  Gen¬ 
eration  2  model  employs  many  of  the  concepts  used 
in  the  Becks tead-Derr-Pr ice  (BDP)  composite  propel¬ 
lant  combustion  model^  with,  however,  major  modi¬ 
fications,  including: 

1.  Variation  in  the  ratio  of  local  oxidizer- 
surface  intersectional  area  to  the  binder  surface 
area  as  the  propellant  surface  regresses  past  an 
oxidizer  particle  is  considered.  (In  the  BDP  model, 
a  geometrical  average  ratio  is  used;  this  involves 
an  assumption  that  a  lot  of  very  nonlinear  processes 
can  be  linearly  averaged.  In  tracking  the  regress¬ 
ion  of  the  flame  front  through  an  oxidizer  crystal 
in  this  first  model  variant,  one  finds  a  strong  var¬ 
iation  in  local  flame  temperature,  stoichiometry, 
and  diffusion  flame  dimensions  between  polar  and 
equatorial  intersections,  as  discussed  in  detail  in 
Ref.  5.) 

2.  The  kinetics  of  subsurface/surface  exo¬ 
thermic  reactions  are  considered,  with  use  of  rate 
expressions  based  upon  the  work  of  Waesche  and 
Wenograd.f?)  (In  the  BDP  model,  subsurface/surface 
heat  release  is  Included  with  the  endothermic  ingred¬ 
ient  vaporization  heats,  with  the  resultant  implicit 
assumption  that  the  amount  of  heat  release  these 
reactions  per  unic  mass  of  propellant  is  independent 
of  such  parameters  as  burning  rate.) 

3.  A  correction  of  an  Inconsistency  in  defi¬ 
nition  of  areas  in  the  BDP  model  is  made. 
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4.  The  calculation  of  the  dimensionless  sto¬ 
ichiometric  group  needed  for  calculation  of  the. 
diffusion  flame  height  via  the  Burke-Schumann^ 
analysis  is  modified.  (The  group  used  in  the  BOP 
model  is  inconsistent  with  that  defined  in  the  orig¬ 
inal  work  of  Burke-Schumann.) 

5.  A  two-flame  (fuel-gas/oxidizer-gas  col¬ 
umnar  diffusion  flame  and  ammonium  perchlorate  mono¬ 
propellant  flame),  rather  than  a  three-flame  model, 
is  used.  (With  correction  of  the  calculation  of 
the  stoichiometry  dimensionless  group  for  the  Burke- 
Schumann  analysis,  it  no  longer  appears  necessary 

to  differentiate  between  the  parts  of  the  diffusion 
flame  inside  and  outside  of  an  anaonium  perchlorate 
monopropellant  flame.)  -(See  Figure  1.) 

6.  The  procedure  for  calculation  of  heat 
feedback  from  the  diffusion  flame  and  the  AP  mono¬ 
propellant  flame  is  modified.  (In  the  BDP  model, 
all  flames  are  considered  to  occur  in  flame  sheets 
at  discrete  distances  from  the  surface:  in  the  cur¬ 
rent  model,  the  AP  monopropellant  heat  release  is 
treated  as  a  flame-sheet  type  heat  release  but  the 
diffusion  flame  heat  release  is  considered  to  occur 
in  a  distributed  fashion  over  a  finite  range  of  dis¬ 
tances  from  the  propellant  surface.)  (See  Fig.  1) 


7.  The  distance  (measured  normal  to  the  pro¬ 
pellant  surface)  associated  with  oxidizer-binder 
gas  interdiffusion  in  the  presence  of  crossflow  is 
assumed  to  be  reduced  by  a  factor,  sin  6,  where  6 
is  the  angle  of  the  resultant  of  the  crossflow  and 
transpiration  velocities  relative  to  the  surface. 
(See  Figure  1.) 


As  mentioned  earlier,  details  of  the  develop¬ 
ment  of  the  first  variant  (baseline)  of  the  Genera¬ 
tion  2  model  are  presented  in  Reference  5.  The 
basic  units  of  the  model  construction  consist  of 
development  of  a  procedure  for  description  of  the 
surface  structure  at  succeeding  incremental  steps 
of  regression  of  the  binder  past  an  oxidizer  par¬ 
ticle,  application  of  an  energy  balance  at  the  sur¬ 
face  along  with  numerous  ancillary  equations  at 
each  increment  for  calculation  of  burning  rate  at 
that  increment,  and  averaging  of  the  incremental 
burning  mass  fluxes  to  obtain  an  overall  average 
burning  rate.  In  this  baseline  variant,  an  organ¬ 
ized  crystal  array  was  assumed  and  used  to  calcu¬ 
late  the  relative  planar  oxidizer  and  fuel  areas  at 
each  intersection  depth.  As  discussed  in  some 
detail  in  Ref.  5,  choice  of  the  best  procedure  for 
calculating  an  average  rate  (the  "end-game"  proced¬ 
ure)  is  not  clear:  several  averaging  equations  are 
presented.  The  prodedure  chosen  for  use  with  the 
first  variant  of  the  model  was: 


avg 


E(A  „  .APOX.  +  mf  .  ,AFU.) 
i  ox,p,j  j  fuel.j  j 

P  _  I  (APOX  -r  AFU.) 

propellant  j  3  j 


(1) 


with  the  summation  being  carried  out  over  all  in¬ 
crements,  including  ones  for  which  the  mass  efflux 
rate  is  zero.  (Since  the  oxidizer  and  fuel  in  gen¬ 
eral  regress  at  different  rates,  it  is  possible  for 
an  oxidizer  particle  to  burn  out  before  the  binder 
reaches  the  bottom  of  its  associated  pit,  leading 
to  aero  burning  rate  for  the  remaining  increments, 
as  discussed  in  Ref.  S.) 


As  explained  in  detail  in  References  2  -  5,  a 
mechanism  by  which  crossflow  is  postulated  by  this 
author  to  affect  the  burning  rate  of  composite 


propellants  involves  the  shortening  of  the  distance 
(measured  normal  to  the  surface)  associated  with  the 
mixing  of  columns  of  fuel  and  oxidizer  gas  leaving 
the  propellant  surface.  In  the  combustion  of  ammon¬ 
ium  perchlorate-loaded  composite  solid  propellants, 
it  is  generally  accepted  that  parallel  columns  of 
oxidizer  and  binder  sublimation/decomposition  pro¬ 
duct  gases  leave  the  surface  from  above  the  oxidizer 
crystals  and  binder,  respectively.  In  the  most  gen¬ 
eral  case  some  heat  is  fed  back  to  the  surface  from 
monopropellant  reaction  of  oxidizer  sublimation  pro¬ 
ducts  while  additional  hear  is  supplied  by  the  mix¬ 
ing  and  reaction  of  the  oxidizer  and  fuel  product 
streams.  (Fig.  1.)  Accordingly,  an  important  fact¬ 
or  in  determining  the  rate  of  heat  feedback  (which 
increases  with  decreased  distance  of  the  gas-phase 
heat  release  zone(s)  from  the  surface)  is  often  the 
rate  of  mixing  of  the  oxidizer  and  binder  gas  pro¬ 
duct  columns.  In  the  absence  of  a  crossflow,  these 
columns  move  perpendicular  to  the  propellant  surface, 
while,  with  crossflow,  they  are  tilted  over  and 
travel  at  an  angle  to  the  surface,  this  angle  being 
determined  by  the  ratio  of  crossflow  velocity  to 
transpiration  velocity  at  any  given  position  above 
the  surface.  As  a  result  of  this  bending,  the  mix¬ 
ing  distance,  measured  perpendicular  to  the  surface, 
decreases.  A  particularly  detailed  description  and 
explanation  of  this  postulated  mechanism  is  presented 
in  Ref.  3  and  4.  In  the  development  of  the  Genera¬ 
tion  1  model  and  in  the  early  stages  of  development 
of  the  Generation  2  model,  this  mechanism  was  assumed 
to  be  the  only  factor  causing  composite  propellant 
burning  rate  to  increase  in  the  presence  of  a  cross- 
flow. 


During  the  past  year,  several  additional  vari¬ 
ants  of  this  second-generation  model  have  been  dev¬ 
eloped  in  a  search  for  an  optimum  model,  both  in 
terms  of  fundamental  physics  and  in  terms  of  explain¬ 
ing  experimental  observations  regarding  the  effects 
of  such  parameters  as  pressure,  oxidizer  particle 
size,  and  propellant  composition  on  burning  rate. 

In  addition,  one  of  the  variants  was  extended  to 
handle  multimodal  oxidizer  formulations.  The  ori¬ 
ginal  erosive  burning  section  of  the  Generation  2 
model  was  found  to  be  inadequate,  as  discussed  later, 
and  this  section  was  subsequently  modified  (in  this 
same  variant)  to  allow  for  augmentation  of  gas- 
phase  transport  properties  by  crossflow-generated 
turbulence  as  a  second  erosive  burning  mechanism. 

The  variants  examined  are  summarized  in  Table  I, 
and  details  of  the  modifications  of  the  original 
baseline  variant  (described  in  detail  in  Ref.  5) 
to  yield  these  latter  variants  are  presented  in 
the  next  section,  along  with  a  description  of  the 
extension  of  the  selected  variant  (IV)  to  handle 
multimodal  oxidizer  cases  and  to  treat  crossflov- 
induced  turbulent  transport  property  augmentation. 
Comparison  of  predictions  and  data  are  then  pre¬ 
sented  in  the  following  section. 

Model  Development 
Unimodal  Oxidizer.  No  Crossflow 

As  indicated  above,  the  first  (baseline)  vari¬ 
ant  of  the  Generation  2  model  is  described  in  detail 
in  Reference  5.  Nomenclature  used  in  the  following 
model  descriptions  is  consistent  with  that  used  in 
that  paper  to  facilitate  comparison.  The  sole  dif¬ 
ference  between  Variant  I  and  Variant  II  lies  in 
replacement  of  the  procedure  of  following  the  reced¬ 
ing  binder  in  incremental  steps  past  the  oxidizer 
crystal  with  use  of  an  "average  burning  surface- 
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oxidizer  intersection"  concept  similar  to  that  used 
by  Beckstead,  Derr,  and  Priced6'  (This  is,  of 
course,  a  major  simplification  to  the  model  since 
the  surface  energy  balance  and  ancillary  equations 
described  in  Ref.  S  need  now  be  solved  only  once 
rather  than  a  number  of  times  equal  to  the  number 
of  incremental  steps  considered  in  the  baseline 
variant,  generally  20  or  more.  Similarly,  the  geo¬ 
metric  relationships  describing  the  surface  only 
have  to  be  solved  once  rather'  than  being  recalcu¬ 
lated  sequentially  at  each  increment.)  As  in  the 
aseline  variant,  the  fuel  is  assumed  to  regress  in 
a  planar  manner  and  the  oxidizer-fuel  surface  is 
forced  to  be  continuous  at  their  intersection. 

These  restrictions,  coupled  with  the  fact  that  the 
linear  regression  rates  of  fuel  and  oxidizer  per¬ 
pendicular  to  their  directions  of  regression  are 
allowed  to  differ,  force  the  oxidizer  to  in  general 
assume  a  curved  shape *as  it  regresses,  thus  creat¬ 
ing  additional  burning  area.  As  in  the  BDT  model, 
upper  hemisphere  and  lower  hemisphere  average  inter¬ 
sections  are  considered  and  values  of  the  curved 
oxidizer  surface  area  (ASOX)  are  calculated  at  each 
intersection  and  averaged  together.  The  "average" 
intersection  in  each  hemisphere  is  defined  such 
that  the  planar  oxidizer-burn  surface  intersectional 
area  (APOX)  is  the  average  of  the  areas  of  all  pos¬ 
sible  intersections  in  that  hemisphere.  With  appli¬ 
cation  of  geometrical  analysis  and  continuity  relat- 


ionships, 

it  is  easily  shown  that: 

APOX  -  itD2/f 
o 

(2) 

AFO  -  (tiDZ/6)(1-VL0)/VL0 

O 

(3) 

In  addition,  the  distances  from  the  initial  top  of 
the  oxidizer  particle  to  the  fuel  surface  at  the 
upper  and  lower  intersections  can  be  shown  to  be 

^Ver  *  <V2)(1  ‘  US** 

(4) 

mT0PI.ower  “  (Do/2)<1  +  lfS*> 

(5) 

Assuming  no  ignition  delay  for  the  oxidizer  parti¬ 
cle,  the  corresponding  distances  from  the  initial 
oxidizer  peak  to  the  center  of  the  oxidizer  crystal 
at  the  times  when  Eqns.  4  and  5  are  satisfied  arc: 

DEL0XUpper“  (rox/rfuel>(XDT0P)Upper  <6) 

^^^Lower  (rox/rfuel)(XD?0P)Lower  (7) 


(with  restraints  that  DELOX  <.  D0) .  Again  through 
geometrical  arguments,  it  may  be  shown  that 


In  Variant  II  of  the  second  generation  model,  Eqns. 
2-9  are  used  to  calculate  AFU  and  ASOX  values 
to  replace  those  used  in  Variant  I.  In  addition, 
the  increment-dependent  flame  temperatures  and  gas 
heat  capacities  of  Eqns.  32  and  33  in  Ref.  S  are 
replaced  by  values  corresponding  to  the  overall 
propellant  stoichiometry.  Finally,  the  burning 
rate  is  calculated  as: 


rox(AS0X) 
ravg  “  APOX 

(replacing  Eqn.  35  of  Ref.  5). 


(10) 


Variant  III  differs  from  Variant  I  in  several 
ways,  all  of  which  are  believed  to  entail  more 
realistic  descriptions  of  various  processes  occur¬ 
ring  during  the  burning  of  a  composite  propellant. 
First,  the  dependence  of  the  O/F  flame  kinetic  dis¬ 
tance,  Lgy,  on  the  flame  temperature  was  changed 
from  that  given  by  Eqn.  25  in  Ref.  5  to: 


RX 
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and  the  corresponding  oxidizer  flame  kinetic  distance 
was  altered  from  Eqn.  26  (Ref.  5)  to: 


AP 


KAP  Vgas.surf  exp  ^EACT,AP7RTAP^ 


.(n-1) 


(12) 


to  allow  dependence  on  local  temperature  (which 
is  calculated  in  a  trial-and-error  loop  on  the  basis 
of  position  of  the  AP  flame  relative  to  the  begin¬ 
ning  and  end  of  the  diffusion  flame).  In  addition, 
the  analysis  of  subsurface  heat  release  was  modified 
to  allow  reaction  only  in  the  liquid  melt  layer  near 
the  surface  and  to  allow  for  reactant  depletion 
(first  order  kinetics)  with  resultant  replacement 
of  Eqn.  20  in  Ref.  5  by: 


a  -  1  -  e-ZINT  (13) 
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(15) 


Finally,  the  assumption  in  Variant  I  that  the  0/F 
flame  heat  release  is  uniform  between  y  =  and 
y  «  +■  FH90sin9  was  replaced  with  a  distributed 

heat  release  more  heavily  weighted  toward  the  base 
of  the  diffusion  flame,  this  weighting  determined 
by  detailed  mixing  profiles  calculated  via  the 

Burke-Schumann Wanalysis.  This  resulted  in  replace¬ 
ment  of  Eqns.  22  and  23  of  Ref.  5  with  considerably 
more  complex  expressions  for  heat  feedback  fluxes 
from  the  gas-phase  flames. 


As  indicated  in  Table  I,  Variant  IV  is  the  same 
as  Variant  III  except  for  use  of  an  "average  inter¬ 
section"  concept:  thus,  modification  of  Variant  III 
to  Variant  IV  follows  sxacfly  the  same  procedure  as 
the  modification  of  Variant  I  to  Variant  II,  des¬ 
cribed  earlier. 


Variants  V  and  VI  represent  modifications  of 
Variant  III  to  allow  different  surface  temperatures 
for  the  binder  and  oxidizer.  To  permit  this,  sepa¬ 
rate  energy  balances  are  written  for  the  oxidizer 
and  binder  (replacing  Eqn.  34  of  Ref.  5),  with  the 
oxidizer  receiving  feedback  from  the  oxidizer  mono¬ 
propellant  flame  and  the  diffusion  flame,  while 
the  binder  receives  feedback  only  from  the  diffusion 
flame.  The  diffusion  flame  feedback  is  apportioned 
between  binder  and  oxidizer  in  the  ratio  of  the 
planar  fuel  and  oxidizer  areas.  Subsurface  heat 
release  is  apportioned  so  that  the  oxidizer  receives 
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■the  fraction,  1/(1  +  SMRBO),  while  the  binder  re¬ 
ceives  the  fraction  SMRB0/(1  +  SMRBO),  where  SMRBO 
is  the  stoichiometric  ratio  (mass)  of  binder  to 
oxidizer.  The  only  difference  between  Variants  V 
and  VI  lies  in  the  treatment  of  the  effect  of  non- 
planar  surface  structure  on  gas  flame  standoff  dis¬ 
tances  relative  to  the  oxidizer  surface.  In  Ver¬ 
sion  V,  allowance  is  made  for  the  non-planar  struct¬ 
ure.  For  oxidizer  protruding,  the  distance  from 
the  oxidizer  surface  to  the  outer  edge  of  the  0/F 
flame  is  reduced  by  the  average  protrusion  height. 

For  oxidizer  recessed,  the  distance  from  the  oxi¬ 
dizer  surface  to  both  the  inner  and  outer  edges  of 
the  0/F  flame  is  increased  by  the  average  recession 
depth.  In  Version  VI,  no  such  adjustment  is  made, 
it  being  assumed  that  the  flame  is  wrinkled  in  such 
a  manner  as  to  track  the  surface  contour. 

Variants  IIIA,  IIIB,  and  IIIC  ail  represent 
modifications  of  Variant  III  aimed  at  increasing 
the  predicted  dependence  of  burning  rate  on  overall 
propellant  flame  temperature.  (As  will  be  discussed 
later,  the  strong  variation  of  flame  temperature 
with  increment  location  in  Variant  III  causes  the 
effect  of  overall  average  flame  temperature  to  be 
fairly  well  washed  out  in  application  of  that  var¬ 
iant.)  Variant  IIIA  is  quire  similar  to  Variant 
III  except  thac  instead  of  using  flame  temperature 
(Tf)  and  gas  heat  capacity  (Cp)  values  corresponding 
to  the  local  oxidizer/fuel  ratio  at  each  increment 
(as  expressed  by  Eqns.  32  and  33  in  Ref.  S),  we  use 
flame  temperature  and  heat  capacity  values  equal  to 
the  overall  propellant  flame  temperature  and  gas 
heat  capacity,  independent  of  increment  number,  at 
each  increment.  That  is,  Eqns.  32  and  33  of  Ref.  5 
are  replaced  by: 

Tf  -  f2  (WFO/ (1-WFO))  (16) 

Cp  «  f3  (WFO/ (1-WFO))  (17) 

where  the  functionalities  are  generated  by  applicat¬ 
ion  of  a  thermochemical  equilibrium  calculation. 
Consistent  with  this  modification,  we  also  slightly 
modify  the  "er.d-game"  procedure  of  calculating  the 
overall  propellant  linear  burning  rate  from  the  mass 
flux  values  calculated  at  each  incremental  inter¬ 
section.  Equation  1  of  tnis  paper  is  still  employed, 
but  the  summations  are  now  limited  to  incremental 
intersections  for  which  me  mass  fluxes  are  non-zero. 
(That  is,  increments  corresponding  to  'ost-oxidizer- 
burnout  are  not  considered,  with  the  assumption 
being  made  that  the  residual  binder  "breaks  off"  in 
some  manner.)  This  has  no  effect  on  the  numerator 
of  Eqn.  1,  but  lowers  the  value  of  the  denominator 
in  cases  where  the  oxidizer  recedes  more  rapidly 
than  the  binder.  (This  modification  of  the  "end¬ 
game"  procedure  is  consistant  with  the  use  of  flame 
temperature  equal  to  overall  propellant  flame  tem¬ 
perature  at  each  incremental  intersection  in  that 
such  use  implies  an  oxidizer/fuel  ratio  in  the  final 
flame  equal  to  the  overall  propellant  oxidizer/fuel 
ratio,  thus  implying  that  the  "left  over"  binder  is 
indeed  somehow  fed  into  the  flame.  It  must  be  em¬ 
phasised  that  this  is  not  a  closed  question,  however. 
As  discussed  in  Ref.  5,  definition  of  an  optimum 
"end-game"  procedure  for  calculation  of  an  average 
linear  propellant  regression  rate  from  values  cal-  . 
culated  at  various  points  in  the  passing  of  the 
burning  front  through  an  oxidizer  particle  is  some¬ 
what  nebulous,  at  best.) 

Variant  IIIB  is  also  a  derivative  of  Variant 
III,  with  arbitrary  assignment  of  fuel  to  oxidizer 
at  various  incremental  intersections  according  to: 


Cx  (D POX) 


-  1  <.  XNA  <  0 


where  XNA  is  an  arbitrary  input  constant  and  is 
determined  for  a  given  value  of  XNA  by  consideration 
of  overall  continuity  for  the  propellant.  A  value 
of  XNA  «  0  corresponds  to  the  ratio  of  fuel  to  oxi¬ 
dizer  pi:  -  intersectional  areas  being  independent 
of  increment  number,  while  XKA"-'0  results  in  polar 
intersections  burning  more  fuel-rich  than  equatorial 
intersections.  As  will  be  discussed  later,  this 
variant  showed  very  little  promise  for  values  of 
XNA  between  0  and  -1,  and  was  quickly  discarded. 

In  Variant"  IIIC,  still  another  offshoot  of 
Variant  III,  the  oxidizer/fuel  ratio  of  gas  streams 
associated  with  a  given  oxidizer  particle  at  each 
ntersection  depth  was  forced  to  be  the  same,  equal 
to  the  overall  propellant  oxidizer/fuel  ratio.  This 
was  accomplished  by  adjusting  the  fuel  area  assoc¬ 
iated  with  the  oxidizer  at  each  intersection  depth 
in  accordance  with: 

(1  -  WFO)  i  „  APOX 

‘"--TBsrsff— 

This  modification  corresponds  to  doing  away  with  the 
concept  of  an  organized  matrix  of  oxidizer  crystals 
in  the  binder  matrix,  employed  in  Variants  I,  III, 

V,  VI,  and  IIIA  (and  described  in  detail  in  Ref.  S), 
and  replacing  it  with  the  assumption  that  random  dis¬ 
tribution  of  oxidizer  crystals  results  in  the  average 
0/F  characteristics  of  a  polar  intersection  (averaged 
over  all  polar  intersections)  and  the  average  0/F 
characteristics  of  an  equatorial  intersection  being 
the  same.  Among  other  things,  this  approach  has 
the  pleasing  effect  of  eliminating  the  problem  of 
"leftover"  binder  which  crops  up  with  the  other  var¬ 
iants.  In  addition,  the  average  linear  propellant 
burn  rates  calculated  using  the  different  averaging 
procedures  described  in  detail  in  Ref.  5  are  more 
nearly  equal  in  this  variant  than  in  Che  other  var¬ 
iants  considering  incremental  passage  of  the  burning 
surface  through  an  oxidizer  particle.  Further,  as 
will  be  discussed  later,  this  variant  is  much  more 
successful  than  Variant  III  in  predicting  dependency 
of  bum  rate  on  overall  propellant  oxidizer/fuel 
ratio  (and  thus  flame  temperature).  A  different 
means  of  calculating  average  burn  rate  from  the 
individual  incremental  rates  (Eqn.  27  in  Ref.  5)  was 
found  to  give  slightly  better  results  than  the  pro¬ 
cedures  used  in  the  other  model  variants.  In  this 
procedure,  the  burning  rate  is  calculated  by  divid¬ 
ing  the  oxidizer  particle  diameter  by  the  sum  of  the 
times  required  for  consumption  of  each  incremental 
thickness  until  the  bottom  of  the  oxidizer  particle 
is  reached :  . 


D  H  TAU, 

0  j  * 


Do/E(A(DELOX)j/rox J 


To  this  author,  this  is  another  plus  point  for  this 
variant  since  Eqn.  20  seems  to  be  more  physically 
realistic  than  Eqn.  1  for  calculation  of  the. average 
linear  regression,  though  it  did  not  produce  satis¬ 
factory  results  in  the  other  variants  studied. 

Multimodal  Oxidizer.  No  Crossflow 

As  will  be  discussed  later.  Variants  IIIA,  IIIC, 
and  IV  were  all  found  to  give  satisfactory  results 
when  tested  against  no-crossflow  burning  rate  data 
for  a  series  of  four  unimodal  oxidizer  composite 
propellants.  Although  Variant  IIIC  is  the  most 
appealing  to  this  author  on  a  physical  basis,  the 
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considerably  lower  computational  complexity  assoc¬ 
iated  with  Variant  IV  led  to  its  selection  for 
extension  for  multimodal  oxidizer  cases.  This  ex¬ 
tension  was  carried  out  in  a  vety  straightforward 
manner  using  Click's  "petit  ensenble"  approach W, 
in  which  a  propellant  containing  oxidizer  particles 
of  different  sizes  is  broken  into  a  series  of  sub¬ 
propellants  or  "pceudopropellants",  each  of  which 
contains  oxidizer  of  only  one  size.  These  sub¬ 
propellants  are  assumed  to  burn  non-interactively, 
with  the  unimodal  oxidizer  model  being  used  to  cal¬ 
culate  a  mass  flux  for  each,  and  straightforward 
averaging  weighted  according  to  fractional  surface 
areas  associated  with  each  subpropellant  then, 
being  used  to  obtain  an  overall  propellant  average 
linear  regression  rate.  The  only  manner  in  which 
oxidizer  of  one  size  is  allowed  to  affect  the  burn¬ 
ing  of  a  subpropellant  containing  oxidizer  of  another 
size  is  through  possible  influence  on  the  assignment 
of  fuel  to  that  subpropellant.  That  is,  rather  than 
fuel  being  assigned  to  each  oxidizer  size  category 
in  direct  proportion  to  the  amount  of  oxidizer  in 
that  category,  the  capability  of  allowing  uneven 
assignment  of  fuel  to  various  oxidizer  size  subpro¬ 
pellants  is  allowed  by  means  of  a  power  law: 


(D)™ 

O  I 


■(21) 


where  Vf is  the  volume  of  fuel  assigned  to  a  par¬ 
ticle  of ’diameter  (D0)it  XEXP  is  an  arbitrary  input 
power  law  constant,  and  C2  is  a  constant  determined 
by  application  of  overall  continuity.  It  may  easily 
be  shown  that  XEXP  »  3  will  result  in  each  subpro¬ 
pellant  having  the  same  oxidizer/fuel  ratio  as  the 
overall  propellant  O/F  ratio.  XEXP<3,  on  the  other 
hand,  will  result  in  subpropellants  with  small  oxi¬ 
dizer  being  more  fuel-rich  than  the  overall  propel¬ 
lant  and  subpropellants  with  large  oxidizer  being 
more  fuel-lean,  with  the  reverse  occurring  for 
XEXP>3. 


Thus,  the  modification  of  Variant  IV  of  the 
Generation  2  burning  rate  model  to  handle  multi¬ 
modal  oxidizer  cases  consists  simply  of  adding  a 
package  at  the  front  of  the  program  to  define  the 
subpropellants,  using  the  existing  progran  to  cal¬ 
culate  burning  rates  for  each  subpropellant,  and 
properly  averaging  these  rates. ^  The  additional 
equations  employed  are  (with  being  defined  as 
the  mass  of  oxidizer  particles  of  particle  diameter 
(D0)i  per  unit  mass  of  total  propellant): 


Modification  of  Erosive  Burning  Package 


Use  of  the  erosive  burning  package  originally 
built  into  the  different  Generation  2  model  variants 
led  in  all  cases,  as  discussed  in  the  next  section, 
to  major  underprediction  of  the  effect  of  crossflow 
on  burning  rate,  Indicating  that  the  proposed  tlame- 
bending  mechanism  was  by  itself  insufficient.  Accord¬ 
ingly,  a  second  possible  mechanism,  augmentation  of 
tu'rbulence  :  ranaport  properties  in  the  region  be¬ 
tween  the  propellant  surface  and  the  gas-phase 
flames  was  invoked  and  combined  with  the  flame-bend¬ 
ing  mechanism.  In  this  approach,  it  was  assumed 
that  both  the  effective  thermal  conductivity  (gov¬ 
erning  feedback  from  the  various  gas  flames)  and 
the  effective  mass  diffusivity  (an  important  para¬ 
meter  in  determining  the  thickness  of  the  diffusion 
flame)  were  increased  in  crossflow  situations  0 y 
crossflow-induced  turbulence.  A  flow  profile  analy¬ 
sis  permitting  calculation  of  eddy  viscosity  (and, 
by  analogy,  total  effective  thermal  conductivity 
and  diffusivity)  as  a  function  of  distance  from 
the  propellant  surface  for  a  given  crossflow  vel¬ 
ocity,  transpiration  velocity  determined  by  the 
propellant  burning  rate) ,  and  temperature  field 
(dependent  on  the  location  of  gas-phase  heat  re¬ 
lease  zones)  was  developed  and  coupled  with  the 
Variant  IV  combustion  model  for  erosive  burning 
calculations.  An  improved  calculation  of  diffusion 
flame-bending  angle  was  also  incorporated  in  this 
analysis. 


For  this  analysis,  preliminary  estimates  of  the 
burning  rate  (and  thus  transpiration  velocity)  and 
the  distance  from  the  surface  to  the  end  of  the  dif¬ 
fusion  flame  (Ljyj  +  FH90  sin  0)  were  first  made 
usvng  laminar  transport  properties  as  in  the  original 
erosive  burning  package.  The  nonblowing  skin  frict¬ 
ion  coefficient  was  then  calculated  as  the  maximum 
of  a  value  calculated  using  a  smooth-wall  equation 
and  one  calculated  using  a  rough-wall  equation: 


cfo  -•  .00140  +  .125  R e-a-**ro.32.  (29) 

Cfo  -  0.95  [4  log10  (R^k)  +  3.48]-2  (30) 


(31) 


ppropellant  1/[Y/P0X+  (l-Y)/ofuell  (23) 

2  •  5  (v<vr“) 

;*»Volume  Fraction  Oxidizer  ,  ... . .XEXP-3, 

1  in  Subpropellant  i  "l/U+eC^D^  105) 

1fJ*Mass  Fraction  Oxidizer  ,  _  ._  ,XEXP-3 , 

1  in  Subpropellant  i  /npoxJ 


was  then  calculated  from  the  previously  estimated 
burning  rate  and  the  freestream  velocity.  This 
parameter  was  next  used  to  estimate  the  actual 
skin  friction  coefficient  with  blowing  (transpirat¬ 
ion).  Unfortunately,  there  is  essentially  no  data 
available  for  values  of  b  greater  than  6-8,  and 
our  range  of  interest  tends  to  go  to  much  higher 
blowing  ratios.  Accordingly,  several  different 
expressions  were  fit  to  existing  data(*0”14)  and 
built  into  the  program  for  optional  means  of  cal¬ 
culating  .skin  friction  with  blowing: 

C./C  =  exp  [-.357  b  +  .0068  b2],  b<  26.25  (32a) 
10  =  ;0092  b>  26.25 

C/C  =  exp  [-.5513b  +  .01316  b2],  b<20  (32b)  . 

1  t0  =  .0031  b>  20 


p|“Density  of 
1  Subpropellant  1 
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(27) 

(28) 


Cf/Cfo  =  exp  (-0.576  b)  (32c) 

Cf/Cfo  =  b/(exp(b)  -  1.0)  (32d) 

Next,  the  wall  shear  stress  was  calculated  as: 
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Application  of  a  momentum  integral  analysis  for  a 
two-dimensional  channel  (the  configuration  in  which 
testing  was  conducted)  yields  the  following  expres¬ 
sion  for  local  shear  stress  as  a  function  of  dis¬ 
tance  from  the  propellant  surface  and  local  cross- 
flow  velocity: 


wall  IW 


K  =  0.9  mTMt  U  S 

INJ  crossflow! 
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'  crossflow 


'  crossflow 


(It  n  -Id  be  noted  that  the  -Ky  term  in  this  equa¬ 
tion  results  from  the  axial  pressure  gradient  tens 
caused  by  friction  and  mass  injection  in  a  confined 
channel:  this  term  has  not  been  included  in  sev¬ 
eral  previous  similar  analyses  of  erosive  burning. 

At  sufficiently  high  blowing  ratios  (leading  to  low 
Twall  values)  this  term  can  be  very  important;  even, 
in  some  cases,  causing  the  calculated  shear  stress 
to  be  negative  away  from  the  surface.  This  author 
interprets  this  condition  as  indicating  boundary 
layer  blowoff:  when  this  condition  is  met,  an 
alternate  boundary  layer  analysis  is  employed,  as 
discussed  later.) 

At  any  rate,  application  of  Eqns.  29-34  permits 
calculation  c  1'  local  shear  stress,  t,  as  a  function 
of  distance  from  the  propellant  surface  (y)  and 
local  crossflow  velocity  (u)  for  a  given  -jet  of  free- 
stream  conditions,  mass  injection  rate,  and  temper¬ 
ature  field.  =  gj  (u,y)l.  An  eddy  viscosity 
approach  is  used  to  relate  the  local  shear  stress 
to  the  local  crossflow  velocity  gradient  via: 

T  =  (u  ‘  pc)  du/dy  (35) 

With  specification  of  an  expression  for  e  as  a 
function  of  y  and  du/dy  (discussed  below) ,  and  use  , 
of: 


P(MW) /RT 


u  =  kT” 


T  ~  Ts  +  I^x+FHSO  sin  0  ",v" 

(the  latter  equation  resulting  from  assumption  of  a 
linear  temperature  profile  from  the  surface  to  the 
end  of  the  diffusion  flame),  Eq.  35  may  be  combined 
with  the  t  -  g^  (u,y)  relationship  obtained  from 
Eqns.  29-34  to  yield  a  first  order  differential 
equation  for  local  crossflow  velocity.  Uith  use  of 
u  =  0  at  y  =  0,  this  is  numerically  integrated  from 
the  wall  to  yield  u(y),  e(y) ,  p(y) ,  and  u(y).  (In 
addition,  with  use  of  v(y)  =  mjNj/o,  v(y)  is  obtain¬ 
ed  and  the  shape  of  streamlines  emanating  from  the 
surface  are  calculated,  permitting  calculation  of 
an  effective  diffusion  flame  bending  angle,  6, 
allowing  for  the  fact  that  the  bent  flame  is  curved 
rather  than  straight,  as  assumed  in  the  earlier 
version  of  the  erosive  burning  package.)  These 
results  are  used  to  calculate 


<Tf  -  yr 


Xeff/XL 


=JZ>eff^£  =  1+0 c/“  *  «2(y' 


That  is,  the  ratio  of  transport  total  properties 
to-  laminar  properties  is  calculated  as  a  function 
of  distance  from  the  surface.-  Average  total  trans¬ 
port  property  values  between  appropriate  zones  are 
then  calculated  and  substituted  for  the  laminar 
values  in  the  diffusional  mixing  equations  and  the 
heat  feedback  equations  in  the  original  model, 
revised  burn  rates  and  flame  distances  are  calcu¬ 
lated,  and  the  procedure  is  repeated  until  converg¬ 


ence  is  achieved.  As  might  be  expected,  this  loop¬ 
ing  procedure  is  considerably  more  complex  in  the 
case  of  multimodal  propellants  than  for  unimodal 
propellants  since  solution  of  the  individual  sub- 
propellant  cases  becomes  interactive  in  the  case 
of  crossflow.  This  interaction  occurs  because  there 
is  only  one  boundary  layer  for  the  overall  propel¬ 
lant  (that  is,  one  cannot  calculate  a  different 
boundary  layer  profile  for  each  subpropellant)  with 
the  boundary  layer  details  being  controlled  by  the 
average  transpiration  velocity,  flame  height,  sur¬ 
face  temperature,  etc.  for  the  overall  propellant 
rather  than  by  the  individual  values  of  these  para¬ 
meters  for  each  subpropellant. 

Several  options  ..or  closure  of  the  boundary 
layer  analysis  through  use  of  an  eddy  viscosity 
equation  were  built  into  the  program:  all  entailed 
use  of  a  Prandtl  mixing-length  type  'expression. 

The  most  comprehensive  of  these  was  a  modified  form 
of  an  emnirical  relation  devr " aped  by  Kays  and 
Moffat which  includes  the  effects  of  blowing 
and  axial  pressure  gradient  l  it  does  not  include 
the  effects  of  roughness.  The  modifications  added 
by  this  author  were  an  attempt  to  include  the  effects 
of  roughness  using  approaches  suggested  by  the  works 
of  van  Driest^®)  and  Cebeci  and  Chang^^/.  The. 
resulting  expression  is:  ' 

-.»W£ 

e=.168 (y+Ay)2  for  ^->65  (40) 

where:  . - 


Ay  =  0  for  - <4.535 


r  =  1+7. 
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The  second  exponential  is  Eqn.  40  is  included  +v 
force  the  damping  function  (major  brackets)  to  go 
to  unity  as  the  roughness  Reynolds  Number  (ku  /v) 
approaches  65  (full  roughness) ,  consistent  with  the 
van  Driest'*"*'  approach.  The  Ay  term  reflects  a  pro¬ 
file  displacement  due  to  roughness  as  suggested  by 
Cebeci  and  Chang(*&) .  Other  options  considered 
included  use  of  the  van  Driest  expression  without 
blowing  or  pressure  gradient  effects,  use  of  the 
above  expression  with  Ay  set  equal  to  zero  for  all 
roughness  values,  and  use  of  an  expression  with  no 
damping  (with  or  without  a  roughness  Ay  factor). 

The  expression  which  thus  far  seems  to  give  the 
best  results  is  the  no-damping,  no-roughness  expres¬ 
sion:  . 

e  =  .168  y*  du/dy  (43) 

though  it  must  be  admitted  that  this  is  not  the  most 
pleasing  on  physical  grounds. 

As  mentioned  earlier,  depending  on  the  Cf/Cf0 
relationship  used  there  are  cases  (high  blowing 
parameter  ones)  where  attempts  to  integrate  the 
shear  stress  equation  out  from  the  surface  results 
in  the  local  shear  stress  going  negative:  this  is 
interpreted  as  representing  some  type  of  boundary 
layer  blowoff.  In  this  case,  it  is  assumed  that 
the  velocity  profile  assumes  the  cosine  law  shape 
(invlscid  no-wall-slip  flow)  measured  by  Yamada 
and  Goto(*7)  among  others,  and  this  profile  is  then 
differentiated  and  the  result  substituted  into  Eqn. 
43  to  yield  the  eddy  viscosity  distribution,  without 
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use  of  momentum  equation,  though  it  is  realized 
that  this  is  not  a  very  satisfactory  resolution  of 
the  problem.  A  major  difficulty  with  the  entire 
boundary  layer  analysis  outlined  above  is  that  the 
bloving  ratio  (transpiration  velocity/crossflow 
velocity)  values  generally  encountered  in  erosive 
burning  are  well  outside  the  range  for  which  bound¬ 
ary  layer  details  have  been  even  marginally  defined. 
To  quote  H.  W.  Liepnann,  "Extrapolation  is  clearly 
anything  but  safe  when  the  methods  used  rely  on 
experimental  confirmation  of  specific  cases  rather 
than  a  deeper  understanding  of  the  fluid  dynamics 
of  turbulence":  such  is  the  case  here. 

Comparison  of  Predictions  and  Data 

A  family  of  six  anroonium  perchlorate  (AP)/ 
hydroxy terminated  polybutadiene  binder  (HTPB)  pro¬ 
pellant  formulations  has  been  used  for  testing  and 
"tuning"  of  the  Generation  2  burning  rate  model 
variants  described  in  the  preceding  section.  Burn¬ 
ing  rate  versus  pressure  data  have  been  obtained 
at  various  crossflow  velocities  (including  zero 
crossflow)  for  these  formulations,  which  are  listed 
in  Table  II.  As  may  be  seen,  four  of  the  formula¬ 
tions  tested  contain  unimodal  oxidizer,  while  the 
other  two  contain  blmodal  oxidizer.  The  unimodal 
formulations  cover  three  particle  sizes  and  two 
oxidizer/fuel  ratios.  In  each  variant  of  the  Gen¬ 
eration  2  model,  there  are  three  "free"  constants 
which  are  optimized  in  terms  of  the  data.  These 
constants  are  Kap» *  pro -exponential  in  the  expres¬ 
sion  relating  the  monopropellant  flame  kinetic 
delay  distance  to  pressure,  temperature,  and  burn¬ 
ing  gas  velocity  (Eq.  12),  KqF,  a  similar  pre-ex¬ 
ponential  for  the  O/F  flame  kinetic  delay  distance 
expression  (Eq.  11),  and  Bsub ,  a  pre-exponential 
for  the  subsurface/surface  reaction  rate  expression. 
The  procedure  which  has  been  employed  in  this  study 
is  to  optimize  these  constants  for  each  variant 
against  the  no-crossflow  bum  rate  data  for  the 
unimodal  propellants,  and  then  leave  them  unchanged 
as  a  given  variant  is  extended  to  multimodal  pro¬ 
pellants  and  crossflow  situations.  Now  that  Variant 
IV  has  been  extended  to  handle  such  cases,  reopti¬ 
mization  of  the  constants  against  all  of  the  data 
for  the  unimodal  and  multimodal  propellants  (with 
and  without  crossflow)  might  be  considered:  at 
this  time,  this  latter  procedure  has  not  been 
employed. 

Unimodal  Oxidizer  Propellants,  No  Crossflow 

Each  of  the  Generation  2  model  variants  was 
tested  against  no-crossflow  data  for  the  unimodal 
propellants  listed  in  Table  II,  with  optimization 
of  the  constants  Kop,  K^p,  and  Bsub.  As  a  result 
of  this  procedure,  several  variants  were  eliminated 
from  further  consideration.  Although  Variants  I 
and  III  appeared  to  provide  comparable  agreement 
with  data.  Variant  I  was  eliminated  because  it  was 
felt  that  the  modifications  made  in  going  from  I 
to  III  resulted  in  a  more  physically  correct  model. 
Similarly,  Variant  II  was  eliminated  relative  to 
Variant  IV. 

With  Variants  III,  V,  and  VI,  values  of  the 
three  free  constants  could  be  found  which  would 
give  excellent  agreement  with  the  no-crossflow 
burn  rate  versus  pressure  data  for  the  first  three 
formulations  of  Table  II:  however,  use  of  these 
same  constants  for  the  fourth  propellant  resulted 
in  serious  underprediction  of  the  burning  rate. 

The  cause  of  this  behavior  was  tracked  down  to  be 


a  tendency  of  these  models  to  underemphasize  the 
effect  of  overall  average  propellant  flame  temper¬ 
ature  on  the  burning  rate.  This  tendency  resulted 
from  the  fact  that  the  averaging  procedure  (for 
averaging  over  the  various  intersection  increments) 
put  a  good  deal  of  weight  on  the  burning  rates  of 
increments  for  which  the  oxidizer/ fuel  ratio  was 
considerably  higher  than  the  overall  oxidizer/fuel 
ratio  (near-equatorial  intersections) .  Since  the 
Vate  of  flame  temperature  increase  with  increasing 
0/F  ratio  decreases  as  the  0/F  ratio  approaches 
stoichiometric,  (all  of  these  formulations  are  fuel- 
rich  overall)  the  flame  temperature  of  the  near- 
equatorial  Intersection  increments  is  nearly  the 
same  for  the  1667®K  formulations  (4525,  4685,  and 
5051)  as  for  the  2065°K  formulation  (5542) .  As  a 
result  the  predicted  burning  rates  for  these  near- 
equatorial  increments  are  nearly  the  same  for  the 
two  20  micron  oxidizer  formulations,  with  the  result 
that  a  set  of  constants  which  gives  good  agreement 
between  theory  and  data  for  the  1667°K  formulations 
tends  to  result  in  underprediction  of  burn  rate 
for  the  hotter  formulation.  It  was  this  result 
that  led  to  development  of  Variants  IIIA,  IIIB,  and 
1IIC,  in  which  more  emphasis  was  placed  on  insuring 
that  overall  formulation  flame  temperature  would 
have  more  direct  influence  on  the  predicted  burning 
rate. 

Variant  IIIB  was  found  to  give  unsatisfactory 
results  for  all  values  of  the  fuel-assignment'  para¬ 
meter  (XNA)  examined  and  was  consequently  dropped. 
The  search  for  an  optimum  value  of  XNA  was  limited 
to  values  between  0  and  -1:  It  may  be  that  a  more 
negative  value  would  result  in  improvement  of  agree¬ 
ment'  between  theory  and  data,  but  it  seemed  more 
worthwhile  to  pursue  Variants  IIIA  and  IIIC  since 
they  seemed  more  attractive  on  physical  grounds 
(particularly  IIIC)  and  required  fitting  of  one 
less  empirical  constant. 

Predicted  and  experimental  burn  rate  versus 
*  'pressure  results  for  the  four  unimodal  formulations 
in  the  absence  of  crossflow  are  presented  for  Var¬ 
iants  IIIA,  IIIC,  and  IV  in  Figures  2,  3,  and  4, 
respectively.  As  may  be  seen,  with  optimization  of 
the  three  free  constants  (slightly  different  for 
the  three  variants,  as  might  be  expected)  all  three 
variants  permit  excellent  prediction  of  the  burning 
rate  versus  pressure  characteristics  of  the  four 
formulations ,  Variant  IIIA  giving  the  best  agree¬ 
ment  of  theory  with  data.  Since  Variant  IV  did  not 
give  appreciably  worse  results  than  the  other  two 
variants  and  was  considerably  simpler  to  extend 
to  treatment  of  multimodal  oxidizer  formulations 
and  to  crossflow  cases  using  the  erosive  burning 
package  described  in  the  previous  section,  it  was 
selected  for  first  extension  to  those  cases. 

Multimodal  Oxidizer  Propellants.  No  Crossflow 

The  extension  of  Variant  IV  of  the  Generation 
2  model  to  handle  multimodal  oxidizer  cases  was 
described  in  the  preceding  section.  Results  of  use 
of  this  extended  model  to  predict  no-crossflow  bum 
r-'.te  versus  pressure  behavior  for  the  two  bimodal 
oxidizer  formulations  listed  in  Table  II  are  pre¬ 
sented  in  Figure  5.  It  should  be  noted  that  the 
values  of  K^p,  Kgp,  and  Bsub  selected  during  analy¬ 
sis  of  the  unimodal  oxidizer  formulations  are  not 
changed:  that  is,  they  are  no  longer  free  constants. 
The  predictions  shown  in  Figure  5  were  made  with 
the  fuel-assignment  constant  (XEXP)  set  equal  to  3 
(resulting  in  equal  stoichiometry  subpropellants). 


As  nay  be  seen,  agreement  between  theory  and  data 
is  excellent  except  at  the  high  pressure  end  of  data 
for  the  90/200  micron  formulation  <5565).  The  two 
highest  pressure  data  points  for  this  formulation 
appear  to  roll  off  severely,  while  the  predicted 
rates  continue  to  climb.  The  reason  for  this  be* 
havior  is  not  clear.  Limited  calculations  using 
lower  values  of  XEXP  (corresponding  to  the  small 
oxidizer  size  subpropollap.t  being  more  fuel-rich 
than  the  large  oxidizer  size  subpropellant)  yielded 
worse  agreement  with  data. 

Erosive  Burning  Predictions  Using  Only  Flame-Bending 

As  discussed  earlier,  all  of  the  model  variants 
examined  originally  contained  an  erosive  burning 
package  based  on  the  assumption  that  the  sole  mech¬ 
anism  causing  burning  rate  increase  with  crossflow 
is  the  bending  over  of  columnar  diffusion  flames  ay 
that  flow.  (This  package  is  described  in  detail 
for  Variant  I  in  Ref.  5.)  For  all  variants,  with 
free  constants  optimized  from  matching  no-cross£low 
data  as  described  earlier,  this  assumption  led  to 
the  underprediction  of  crossflow  effects  for  all 
formulations  at  all  pressures  and  crossflow  veloci¬ 
ties.  An  example  of  this  underprediction  is  shown 
in  Figure  6,  where  predictions  made  using  Variant 
IV  with  the  original  erosive  burning  package  are 
compared  with  data  for  Formulation  4525  (73/27  AP/ 
HTPB,  20  micron  AP).  Based  on  these  results,  the 
erosive  burning  portion  of  the  Generation  2  model 
was  revised  as  described  in  the  previous  section  to 
allo<  consideration  of  augmentation  of  gas  trans¬ 
port  properties  by  crossflow-induced  turbulence  as 
well  as  flame  bending. 

Erosive  Burning  Data  Comparison  With  Revised  Model 
Predictions 

Several  different  eddy  viscosity  models  and 
several  equations  for  extrapolation  of  skin  friction 
data  to  high  blowing  ratio  conditions  were  built 
into  the  revised  Variant  IV  model,  as  discussed 
earlier.  Studies  to  date  indicate  that  the  erosive 
burning  data  for  the  formulations  listed  in  Table  II 
are  best  predicted  using  the  no-damping,  no-rough¬ 
ness  eddy  viscosity  expression  (Eq.  43)  and  the 
first  of  the  four  Cf/Cf0  expressions  (Eq.  32a),  but 
a  firm  conclusion  regarding  the  most  suitable  ex¬ 
pressions  Has  not  yet  been  made.  (As  indicated, 
a  major  difficulty  with  this  analysis  is  that  it 
involves  making  boundary  layer  predictions  in  a 
blowing-ratio  regime  which  has  not  been  studied, 
and  extrapolation  of  turbulent  flow  field  d3ta  is 
extremely  risky.)  Theoretical  predictions  and  data 
for  the  six  formulations  tested  to  date  are  pre¬ 
sented  in  Figures  7-12.  As  may  be  seen,  the 
agreement  between  theory  and  data  for  five  of  the 
six  formulations  is  quite  encouraging.  With  Formu¬ 
lation  5565  (the  82/18  AP/HTPB  bimodal  formulation 
with  90  and  200  roirron  ammonium  perchlorate)  how¬ 
ever,  the  model  predicts  considerably  more  sensi¬ 
tivity  of  burning  rate  to  crossflow  than  actually 
observed.  Use  of  Eqn.  32d  for  calculation  of  Cj/ 

Cf0  results  in  considerable  improvement  of  agree¬ 
ment  between  theory  and  data  for  this  formulation, 
but  only  at  the  expense  of  worsening  agreement 
between  theory  and  data  for  the  other  five  formu-  ' 
la t ions.  Thus  far,  no  combination  of  eddy  viscos¬ 
ity  equations  and  Cf/C f0  expressions  considered  has 
resulted  in  improvement  of  the  predictions  for  the 
"problem"  formulation  without  simultaneously  de¬ 
grading  the  agreement  between  theory  and  experiment 
for  the  other  formulations. 


Considerable  thought  has  been  devoted  to  this 
problem,  resulting  in  definition  of  an  interesting 
dilemma.  Formulation  4525  (73/27  AP/HTPB,  1667'K 


flame  temperature,  20  micron  oxidizer)  and  Formu¬ 
lation  5565  (82/18  AP/HTPB,  2575°K  flame  temperature 
bimodal  mix  of  90  and  200  micron  oxidizer)  were 
tailored  to  give  nearly  identical  burn  rate  versus 
pressure  characteristics  in  the  absence  of  cross- 
flow.  This  near  equality  of  burning  rates  resulted 
from  cancellation  of  a  considerably  higher  temper¬ 
ature  driving  force  (Tfiame-  Tsur£acc)  for  the 
latter  formulation  by  a  lower  diffusion  flame  dis¬ 
tance  associated  with  the  smaller  oxidizer  formula¬ 
tion.  That  is,  in  a  no-crossflow  situation  the 
diffusion  flame  zone  stretches  considerably  further 
from  the  surface  for  the  5565  formulation,  offset¬ 
ting  the  effect  of  higher  flame  temperature  on  burn 
rate.  Thus,  it  would  appear  that  the  heat  release 
zones,  being  much  further  out  into  the  higher  tur¬ 
bulence  (and  stronger  flame-bending)  region  of  the 
boundary  layer  for  the  5565  formulation  should  be 
much  more  affected  by  crossflow  for  this  formula¬ 
tion  than  for  Formulation  4525,  and  thus  the  burn¬ 
ing  rate  of  5565  should  be  much  more  sensitive  to 
crossflow  than  that  of  4525.  Indeed,  this  is  what 
is  predicted  by  the  model  and  is  what  it  seems  to 
this  author  must  be  predicted  by  any  model  based 
on  either  diffusion- flame  bending  or  turbulence 
augmentation  of  transport  properties.  However, 
comparison  of  the  data  in  Figures  8  and  12  reveals 
that  the  erosive  burning  characteristics  of  these 
two  formulations  are  essentially  identical.  This 
appears  to  be  a  fundamental  problem  for  which  this 
author  does  not  currently  have  an  explanation. 


Summary 

Several  variants  of  a  model  for  the  non-erosive 
burning  of  unimodal  oxidizer  composite  propellants 
have  been  developed.  These  variants  have  been 
tested,  with  optimization  of  three  "free"  constants 
against  no-crossflow  burning  rate  data  for  a  family 
of  four  AP/HTPB  propellants  experimentally  studied 
at  Atlantic  Research.  Three  variants,  including 
two  which  treat  the  increment  -1  regression  of  the 
burning  front  through  an  oxidizer  crystal  and  one 
in  which  an  average-intersection  concept  is  employed 
have  been  found  to  predict  the  data  quite  well. 

The  average-intersection  variant  was  extended  to 
handle  crossflows  and  multimodal  oxidizer  formula¬ 
tions.  This  extended  model  was  found  to  give  excel¬ 
lent  agreement  with  no  crossflow  burning  rate  data 
for  two  additional  formulations  containing  bimodal 
oxidizer.  In  the  initial  development  of  the  erosive 
burning  aspect  of  the  model,  columnar  diffusion 
flame  bending  was  assumed  to  be  the  only  mechanism 
leading  to  bum  rate  augmentation  by  crossflows. 

This  assumption  led  to  severe  underprediction  of 
erosive  burning  effects,  and  a  flow  profile  analysis 
for  prediction  of  turbulence  augmentation  of  trans¬ 
port  properties  as  well  as  flame-bending  was  sub¬ 
sequently  built  into  the  model.  This  revised  model 
is  found  to  give  good  agreement  with  erosive  burn¬ 
ing  data  for  five  of  the  six  formulations  studied, 
but  gives  predicted  rates  higher  than  observed  for 
the  sixth  formulation.  It  is  felt  that  this  problem 
results  from  an  inability  to  accurately  predict 
turbulence  profiles  at  high  blowing  ratios  (trans¬ 
piration  velocity/crossflow  velocity)  in  combinat¬ 
ion  with  the  fact  that  the  formulation  in  question 
contains  very  large  oxidizer  particles,  making  its 
predicted  burning  rate  quite  sensitive  to  the 
boundary  layer  details < 
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Nomenclature 


tuel  surface  area  associated  with  APOX 
planar  projection  ot  exposed  oxidizer 
particle  surface  area 
total  curved  oxidizer  exposed  surface 


eACT,OF 


eACT,AP 


lKX 

Crossflow 

MW 

®INJ 


Ucr oss flow 

Wfs 

u* 

«L 


temperature 

temperature  at  location  of  monopropel¬ 
lant  oxidizer  flame 
flame  temperature 
oxidizer  melting  temperature 
propellant  bulk  (initial)  temperature 
surface  temperature 

time  for  fuel  regression  plane  to  move 
from  one  increment  to  the  next 
crossflow  velocity  at  distance  y  from 
the  surface 

mean  crossflow  velocity 
mainstream  crossflow  velocity 
friction  velocity,  Ufs  \fctl2 
defined  by  Eq.  IS 


vgas,surf 


pre-experimental  for  subsurface  reaction 
rate  equation 

blowing  parameter,  defined  by  Eq.  31. 
skin  friction  coefficient,  including 
transpiration  effects 
no-blowing  skin  friction  coefficient 
gas  heat  capacity 
solid  oxidizer  heat  capacity 
oxidizer  particle  diameter 
distance  of  center  of  oxidizer  crystal 
surface  from  initial  oxidizer  peak 
diameter  corresponding  to  APOX 
gas  diffusivity  (molecular) 
effective  total  gas  diffusivity  (mole¬ 
cular  plus  turbulent) 
activation  energy  for  fuel-oxidizer 
gas  reaction 

acitvation  energy  for  monopropellant 
oxidizer  gas  reaction 

activation  energy  for  subsurface  reaction 

distance  required  for  mixing  of  oxidizer 

and  fuel  gas  screams 

channel  half-height 

roughness  height 

constant  defined  by  Eq.  34 

constant  in  expression  for  oxidizer 

monopropellant  reaction  distance 

constant  in  expression  for  0/F  gas 

reaction  distance 

oxidizer  monopropellant  gas  reaction 
distance 

oxidizer/fuel  gas  reaction  distance 
cross  flow  Mach  number  (mean) 
molecular  weight 

surface  blowing  mass  flux  (burning 
mass  flux  based  on  planar  area) 
mass  flux  of  fuel 

oxidizer  mass  flux,  based  on  planar 
surface  projection 

mass  flux  of  propellant  (or  subpro'pel- 
lant)  based  on  planar  surface  area 
global  gas-phase  reaction  order 
pressure 

gas  law  constant 

Reynolds  Number  based  on  channel  hydrau¬ 
lic  diameter 

average  linear  regression  rate  of  pro¬ 
pellant 

linear  regression  rate  of  fuel  surface 
linear  regression  rate  of  oxidizer, 
normal  to  its  curved  surface 


transpiration  velocity  at  distance  v 
from  the  surface 

volume  of  fuel  associated  with  an  oxidizer 

particle  of  diameter  di 

gas  velocity  away  from  the  propellant 

surface 

volumetric  oxidizer  fraction  in  propellant 
weight  fraction  oxidizer  in  propellant 
distance  of  fuel  surface  from  initial 
oxidizer  peak 

exponent  used  is  assignment  of  fuel  to 
various  oxidizer  particle  size  classes 
exponent  used  in  assignment  of  fuel  sur¬ 
face  area  to  various  oxidizer-surface 
intersections  (Eq.  18) 
distance  from  propellant  surface 
defined  by  Eq.  14 

mass  fraction  of  oxidizer  reacting  at 
or  below  propellant  surface 
mass  fraction  oxidizer  in  overall  propel¬ 
lant 

mass  of  oxidizer  of  particle  size  class  1 
per  mass  of  propellant 

mass  fraction  oxidizer  in  subpropellant  i 
specific  heat  ratio 

offset  distance  due  to  roughness  (Eq.  41) 
eddy  viscosity 

volume  fraction  oxidizer  in  subpropel¬ 
lant  i 

effective  total  gas  thermal  conduct¬ 
ivity  (molecular  plus  turbulent) 
molecular  gas  thermal  conductivity 
condensed-phase  oxidizer  thermal  con¬ 
ductivity 

(molar  fuel/oxidizer  ratio) /stoichio¬ 
metric  fuel/oxidizer  ratio) 
flow-bending  angle  (See  Fig.  1) 
density  of  subpropellant  i 
oxidizer  density 


Ppropellitit  propellant  density 


fuel  (binder)  density 
gas  density 

mainstream  gas  density 

viscosity  (laminar) 

shear  stress  at  distance  y  from  the 

propellant  surface 

wall  shear  stress 

kinematic  viscosity,  p/p 
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Table  I.  Unimodal  Oxidizer  (No  Cross  Flow)  Generation  2  Model  Variants  Examined 

I.  ORIGINAL  GENERATION  2  MODEL.  DESCRIBED  IN  DETAIL  IN  REFERENCE  5.  UNIFORM 
DISTRIBUTION  OF  O/F  FLAME  HEAT  RELEASE  BETWEEN  X  »  Lrx  AND  X  •  Lrx  ♦  FH90.  LAP 
INDEPENDENT  OF  LOCAL  TEMPERATURE.  STEP  THROUGH  VARIOUS  INTERSECTIONS  OF  OXIDIZER 
AND  FUEL.  WITH  3  WAYS  OF  AVERAGING  BURN  RATE  EMPLOYED. 

II.  SAME  AS  I.  EXCEPT  FOR  USE  OF  AN  AVERAGE  FUEL-OXIDIZER  INTERSECTION  CONCEPT. 

III.  MODIFICATION  OF  I.  ALLOWANCE  FOR  DEPENDENCE  OF  LAP  ON  LOCAL  TEMPERATURE.  ALTERED 
DEPENDENCY  OF  LRX  ON  Tf|ame-  SUBSURFACE  HEAT  RELEASE  CALCULATION  MODIFIED  TO  ALLOW 
FOR  REACTANT  DEPLETION.  WEIGHTING  OF  DISTRIBUTION  OF  O/F  FLAME  HEAT  RELEASE  BETWEEN 
TOP  ANO  BOTTOM  OF  CONICAL  FLAME  CHANGED  TO  BE  MORE  CONSISTENT  WITH  DETAILED  MIXING 
PROFILES  CALCULATED  VIA  BURKE-SCHUMANN  ANALYSIS. 

IV.  SAME  AS  III.  EXCEPT  FOR  USE  OF  AN  AVERAGE  FUEL-OXIDIZER  INTERSECTION  CONCEPT. 

V.  SAME  AS  III,  EXCEPT  FOR  ALLOWANCE  OF  DIFFERENT  SURFACE  TEMPERATURES  FOR  BINDER  AND 
OXIDIZER.  DIFFUSION  FLAME  HEIGHTS  RELATIVE  TO  AP  SURFACE  ADJUSTED  TO  ALLOW  FOR  ITS 
PROJECTION  OR  DEPRESSION  RELATIVE  TO  BINDER  CONTINUUM  SURFACE. 

VI.  SAME  AS  V.  EXCEPT  NO  ADJUSTMENT  OF  DIFFUSION  FLAME  HEIGHTS  {WRINKLED  FLAME). 

1 1 1  A.  SAME  AS  III.  EXCEPT  THAT  FLAME  TEMPERATURE  IS  SET  EQUAL  TO  FLAME  TEMPERATURE 

ASSOCIATED  WITH  OVERALL  COMPOSITION  AT  ALL  INTERSECTIONS. 

IIIB.  SAME  AS  III.  EXCEPT  ASSIGNMENT  OF  FUEL  AREA  TO  EACH  INCREMENTAL  INTERSECTION  IS  DONE 
VIA  AFU/APOX  *  C(00)XNA. 

Hit  SAME  AS  III.  EXCEPT  O/F  RATIO  IS  FORCED  TO  BE  CONSTANT  (AT  OVERALL  PROPELLANT  O/F)  AT 

EACH  INTERSECTION  BY  ADJUSTMENT  OF  THE  ASSOCIATED  FUEL  AREA  AS  THE  OXIDIZER  - 
SURFACE  INTERSECTIONAL  AREA  CHANGES. 


Table  It.  Propellants  Used  in  Model  Checkout. 
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FORMULATION 
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DESIGNATION 
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(a)  No  Cross  Flow  Case  (b)  Cross  Flow  Case 

Figure  t.  Schematic  of  Burning  Composite  Propellant,  With  and  Without  Cross  Flow 


BURN  RATE  tin/Mc)  BURN  RATE  (in/tac) 


■■■■I 

■■■II 

Hill 

■ini 


mmuzkw 

mtp&Mw*01 

wzmm 


0 - 

—  5/JAP.  73/27 

o - 

—  20/JAP.  73/27 

A - 

—  200/JAP.  73/27 

i  t 

—  20/JAP.  77/23 

1  1  1  1  II 

20  40  CO  M  100 

PRESSURE  (xml 


msimuim 


BSUB  •  8.70E  12 


40  60  «0  100 

PRESSURE  (atm) 


Figure  2.  Model  Variant  IIIA  -  Theory  and  Oata  for  Four ' 
Unimodal  Oxidizer  Formulations 


Figure  3.  Model  Variant  IIIC-Theory  and  Data  for  Four 
Unimodal  Oxidizer  Formulations 
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Figure  4.  Model  Variant  IV  -  Theory  and  Oata  for  Four 
Unimodal  Oxidizer  Formulations 


Figure  5.  Model  Variant  IV  -  Theory  and  Oata  for  Two 
Multimodal  Formulations 


ftATC  (tm/ttc 


BURNING  RATE  <cm/Mc) 


3B aaBSBEEI.  - ■  ■'  --  - 


Figure  1 0.  Erosive  Burning  Data  and  Predictions  -  77/23 
AP/HTPB,  20jrAP 


Figure  11.  Erosive  Burning  Data  and  Predictions  -  82/18 
AP/HTPB.  41%  IfiAP,  41%  7/iAP 


Figure  12.  Erosive  Burning  Date  and  Predictions  -  82/18 
AP/HTPB.  13.65%  90pAP,  68.35%  200jrAP 


